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ABSTRACT 

As massive black holes (MBHs) grow from lower-mass seeds, it is natural to expect 
that a leftover population of progenitor MBHs should also exist in the present universe. 
Dwarf galaxies undergo a quiet merger history, and as a result, we expect that dwarfs 
observed in the local Universe retain some 'memory' of the original seed mass distri- 
bution. Consequently, the properties of MBHs in nearby dwarf galaxies may provide 
clean indicators of the efficiency of MBH formation. In order to examine the properties 
of MBHs in dwarf galaxies, we evolve different MBH populations within a Milky Way 
halo from high-redshift to today. We consider two plausible MBH formation mech- 
anisms: 'massive seeds' formed via gas-dynamical instabilities and a Population HI 
remnant seed model. 'Massive seeds' have larger masses than PopIII remnants, but 
form in rarer hosts. We dynamically evolve all halos merging with the central system, 
taking into consideration how the interaction modifies the satellites, stripping their 
outer mass layers. We compare the population of satellites to the results of N-body 
simulations and to the observed population of dwarf galaxies. We find good agreement 
for the velocity, radius and luminosity distributions. We compute different properties 
of the MBH population hosted in these satellites. We find that some MBHs have been 
completely stripped of their surrounding dark matter halo, leaving them 'naked.' We 
find that for the most part MBHs retain the original mass, thus providing a clear 
indication of what the properties of the seeds were. We derive the black hole occupa- 
tion fraction (BHOF) of the satellite population at 2 = 0. MBHs generated as 'massive 
seeds' have large masses that would favour their identification, but their typical BHOF 
is always below 40 per cent and decreases to < 1 per cent for observed dwarf galaxy 
sizes. In contrast, Population HI remnants have a higher BHOF, but their masses have 
not grown much since formation, inhibiting their detection. 

Key words: galaxies: dwarf - galaxies: evolution - galaxies: formation - black hole 
physics - cosmology: theory 



1 INTRODUCTION 

Deriving observational constraints on MBH formation, oc- 
curring at very high redshifts, is clearly challenging. The 
first 'seeds' could be as light as a few hundred Mq. Also, 
initial conditions tend to be erased very fast if accretion is ef- 
ficient in growing MBHs (^Volontcri & Gnedin 2009). Galac- 
tic archeology offers an alternative to high-redshift observa- 
tions. Simple arguments suggest that MBHs might inhabit 
also the nuclei of dwarf galaxies, such as the satellites of the 
Milky Way and Andromeda, today. As MBHs grow from 
lower-mass seeds, it is natural to expect that a leftover pop- 
ulation of progenitor MBHs should also exist in the present 
universe. Indeed, we expect that one of the best diagnostics 
of 'seed' formation mechanisms would be to measure the 
masses of MBHs in dwarf galaxies. The progenitors of mas- 



sive galaxies have in fact a high probability that their central 
MBH is not 'pristine', that is, it has increased its mass by 
accretion, or it has experienced mergers and dynamical in- 
teractions. Any dependence of the MBH mass, Mbh, on the 
initial seed mass is largely erased. In contrast, dwarf galax- 
ies undergo a quieter merger history, and as a result, at low 
masses the MBH occupation fraction and the distribution 
of MBH masses still are expected to retain some 'memory' 
of the original seed mass distribution. The signature of the 
efficiency of the formation of MBH seeds will cons equently 
be stronger in dwarf galaxies (IVolonteri et al.ll2008l ). 

The records for the smallest known MBHs belong to 
the dwarf Seyfert 1 galaxies POX 52 and NGC 4395. They 
are believed to contain black holes of mass Mbh ~ 10'' M© 
(|Barth et al.ll2004l : [Petersen et al.ll2005h . There are also sig- 
nificant non-detections of MBHs in a few nearby galaxies 
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from stellar-dynamical observations, most notably the Local 
Group Scd-type spiral galaxy M33, in whic h the upper limit 
to MBH is just a few thou sand solar masses ijCebhardt et alJ 
I2OOII : iMerritt et all I2OOII '). Similarly, in the Local Group 
dwarf eUiptical gal axy NGC 205, A/bh < 3.8 x 10* M© 
l|Valluri et al.ll2005h . These results suggest that the MBH 
'occupation fraction' in low-mass galaxies might be signifi- 
cantly below unity, but at present it is not possible to carry 
out measurements of similar sensitivity for galaxies much 
beyond the lim i ts of the Local Group (llbata et al 1 12OO9I : 
iLora et al ]|2009l : iGreene fc h"c1|2007I : iDebattista et al.ll2006h . 

In this paper we explore theoretical expectations for (1) 
the probability that a dwarf galaxy hosts an MBH, and (2) 
if the merging and accretion history in dwarf galaxies leads 
to different scaling relationships of MBHs with their hosts. 
We explore here the predictions of different models of MBH 
'seed' formation, and derive constraints that may be testable 
with current and future instruments. 



2 MASSIVE BLACK HOLE FORMATION AND 
DYNAMICAL EVOLUTION 

We follow the formation and evolution of an MBH pop- 
ulatio n in a Milky Way -size halo in a ACDM uni- 
verse ijDunklev et al. I2OO8I ). Our technique follows that of 
IVolonteri et all (|2003h . as we use Monte Carlo realisations 
of the merger histories of dark matter halos. We analyse 
here 5 different realisations of halos that reach a mass of 
Mh = 2 X 10^^ Mq at 2 = 0. We seed the high-redshift 
progenitor halos with black holes and follow them from for- 
mation to 2 = 0. We focus on the signatures of black hole 
formation efficiency in satellite galaxies surviving until to- 
day. 



2.1 Massive black hole formation models 

We adopt three different models to seed haloes with black 
holes: massive seeds and two Population HI remnant seed 
models. These models determine which haloes are seeded 
with black holes at high redshifts. They also set a minimum 
mass for central MBHs today. 



2.1.1 Population in remnants 

For the Population HI remnant models, we follow two 
schemes that differ only in the efficiency of MBH seed for- 
mation. We assume here that MBHs form as end-product 
of the very first generation of stars. The main features of 
a scenario for the hierarchical assembly of MBHs left over 
by the fir st stars in a ACDM c osmology have be e n dis - 
cussed by 'Volonteri et al.' ('200^ IVolonteri fc ReesI (|2006l ') 
and IVolon tcri fc Natarajan (.20091 '). The first stars are be- 
lieved to form at 2 ~ 20 — 30 in minihalos, Afmin ~ 10® M©, 
above the cosmological Jeans mass. They collapse at 2 ~ 
20 — 50 from the rarest v-a peaks of the primordial den- 
sity field. In this regime, cooling is possib le by means of 
molecular hydrogen (Tyir > 2 — 3 x 10^ K, iTegmark et al.l 
ll997l : lYoshida et al.|[200^ , but the inefficient cooling at zero 
metallicity might lead to a very top-heavy initial stellar mass 
function. Specifically, the e arliest-forming s t ars are likely 
to have been very massive (jCarr et al ] 1 1984 iBromm et al.l 



1 19991 , 12OOII : lAbel et alllioool : lYoshida et aUbOOSl ). If stars 
form above 260 Mq , they woul d rapidly collapse to massive 
black holes with little mass loss (iFrver et aLlbOOll ). i.e., leav- 
ing behind seed MBHs with masses Mbh ~ 10^ - 10^ Mq 
(Madau & Rees 2001). We here consider that MBH seeds 
populate haloes with formation redshift 2 > 12 which rep- 
resent density peaks with i^c = 3 ('peak3'; as a reference: 
Afh > lO'^ Mq at 2 = 20, and Mh > 10* Mq at 2 = 12) 
or Vc = 3.5 ('peaks. 5'; as a reference: Afh > 10® Mq at 
2 — 20, and Afh > 10^ Mq at 2 = 12), while also requiring 
that Tvir ^ 2500 K to ensure effective H2 cooling (imply- 
ing masses above 3 x 10® Mq). Seeds form as 100 Mq black 
holes. 



2.1.2 Massive MBH seeds 

The massive seed model relies instead on the col- 
lapse of supermassive objec ts for med directly o ut of 

dense gas (Koushiappa s et al.l |2004| : iBegelman et all I2OO6I : 
[Eodato & Natarajan 

2ooir 

and references therein), where 
the mass infiow is regulated by the degree of stability of 
the cooling gas. Here we assume that gas is accumulated 
in the centre of a halo via viscous instabilities. A dynam- 
ically unstable disc can develop non-axisymmetric spiral 
structures that effectively redistribute angular momentum, 
causing mass inflow. This process stops when the amount 
of mass transported to the centre is enough to make the 
disc marginally stable. We here follow iLodato fc NataraianI 
(200^), who suggest that the mass inflow can be computed 
from the Toomre stability criterion and from the disc proper- 
ties, determined from the dark matter halo properties (halo 
mass, Mh, virial temperature, Tvir an d spin parameter. A; 
I Mo et al. 199 8). We ref e r the reader to lLodato fc NataraianI 
l|2006l ): IVolonteri et al.l (Hool) for a comprehensive descrip- 
tion. We summarize here the main features of the model, 
and how we operationally implement it. 

Consider a dark matter halo of mass Afh and virial tem- 
perature Tvi£], containing a gas mass in cold gas Afgas = 
fdMh (we assume that th e gas fraction co oling in is roughly 
5% implying /d = 0.05, IMo et al.|[l993 ). The other main 
parameter characterizing a dark matter halo that is rele- 



vant here is its spin parameter A (= JhE^^^^ /GM^^^ , where 
Jh is the total angular momentum and Eh is the binding 
energy). The distribution of spin parameters for dark mat- 
ter halos measured in numerical simulations is well fit by 
a lognormal distribution in Aspin, with mean Aspin = 0.05 
and s tandard deviation tJx = 0.5 (e.g.. [ Bullock et al.ll200ll : 
Ivan d en Bosch ot al. 20o3; iMaccio et aLlbOOSl ). 

If the virial temperature of the halo Tvir > 2gas, the gas 
collapses and forms a rotationally supported disc. For metal- 
free gas, the cooling function is dominated by hydrogen. In 



f5/2 



^ A halo at redshift z is uniquely characterized by a virial radius 
rvir, defined as the radius of the sphere encompassing a mean mass 
overdensity (5vir- ^vir ~ 178 in an Einstein-de Sitter Universe. De- 
tailed cal culations for d i fferent cos mologies (e. g ., AC DM) can be 
found in iLacev fc Cold lll993D and lEke etal] lll996l) . Prom the 
virial theorem the virial mass, Mh can be calculated straightfor- 
wardly, along with the circular velocity, Vc = y/GM\Ji\i^ and 
virial temperature T^ir = ^J-mpV^ / Ciks) , where ^ ~ 0.722 is 
the mean molecular weight, rrip is the proton mass and ks is the 
Boltzmann constant. 
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thermal equilibrium, if the formation of molecular hydrogen 
is su ppressed (see the discussion in iDevecchi fc Volonteril 
I2OO9I ). these discs are expected to be nearly isothermal 
at a temperatu re of a few thousand Kelv in (here we take 
Tgas « 5000K, iLodato fc NataraianllioO^ ). For low values 
of the spin parameter A, the resulting disc can be compact 
and dense and is subject to gravitational instabilities. This 
occurs when the Toomre stability parameter Q approaches 
a critical value Qc of order unity (following Volonteri et al. 
2008 we adopt Qc = 2, in order to match observational con- 
straints on the MBH and quasar population). If the destabi- 
lization of the system is not too violent, instabilities lead to 
mass infall instead of fragmentation in to bound clumps and 
globa l star formation in the entire disc (|Lodato fc NataraianI 
l2006h . This is the case if the inflow rate is below a critical 

threshold Mmax = StJc-jf that the disk i s able to sustain 
(where Ofc ~ 0.12 describes the viscosity. iRice et al.ll2005l ) 
and molecular and metal cooling are not important. 

This mass redistribution process ceases when the 
amount of mass transported to the center, Ala, is enough 
to make the disc marginally stable. This can be computed 
easily from the Toomre stability criterion and disc proper- 
ties, determi ned from the da rk matter halo mass and angular 
momentum dMo^eTaLlllHgi): 



Ma = <^ 
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T 
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where 

Amax — 



/dQc/8(/d/jd)(rvi./rgas) 



1/2 



A < A 

max 

A > A max 
(2) 



is the maximum halo spin parameter for which the disc is 
gravitationally unstable, and jd is the fraction of the halo 
angular momentum retained by the collapsing gas {jd ~ fd 
if specific angular momentum is conserved). 

For large halo masses, the internal torques needed to 
redistribute the excess baryonic mass become too large to be 
sustained by the disc, which then undergoes fragmentation. 
This occurs when the virial temperature exceeds a critical 
value Tmax, given by: 

2/3 

(3) 



f 4ac 



V /d 1 + Msii/fiMh 



To summarize, every dark matter halo is characterized 
by its mass Mh (or virial temperature Tvir) and by its spin 
parameter A. The latter is drawn from a lognormal distribu- 



tion in Aspin with mean As 



: 0.05 and standard deviation 



ax = 0.5 iMaccio et al.l [20071 . and references therein). The 
gas has a temperature Tgas = 5000K. If A < Amax (see Eq[2]) 
and Tvir < Tmax (Eq. [21), then we assume that a seed BH of 
mass Mbh = Ma given by eqn. IT} forms in the center. 

The gas made available in the central compact re- 
gion can then form a central massive object, for in- 
stance via the intermediate stage of a 'supermassiv e' star 
iHovle fc Fowleil Il963l : iBaumgarte fc Shapiro! 1 19991 ). or a 
'quasistar' (an initially low-mass black hole rapidly accret- 
ing with in a massive, radiation-press ure-supported envelope, 
see also iBegel man et al.ll2006l . l2008h . Hence, the black hole 
seed mass estimates based on Equation [T] should be con- 
sidered as upper limits. We here consider that MBH seeds 
with Mbh — Ma. can form in haloes with formation redshift 



2; > 12 that satisfy all the above criteria. Th e mass function 
of see ds peaks at 10^ — 10^0 (see Figure 2 in I Volonteri et al.l 
l2008l ). 

We remind the reader here that this process is effective 
for halos with low angular momentum (low spin parame- 
ter, A < 0.01) and zero metallicity, where cooling is driven 
by atomic hydrogen cooling and the difference between gas 
and virial temperature is small (making the disc resilient 
to global fragmentation and star formation). The efficiency 
of the seed assembly process ceases at large halo masses 
(Tvir > 1.4 X 10'' K), where the mass-accretion rate from the 
halo is above the critical threshold for fragmentation and 
the disc undergoes global star formation instead. 



2.2 Massive black hole and galaxy evolution 

To study the MBHs at 2 = 0, we follow the evolution of the 
black holes along with their host galaxies using models that 
track both black hole mergers and accretion. Several recent 
works poin t out how fragile the environment of low-mass 
galaxies is l|Bovill fc Ricottil |2009| . and references therein). 
The shallow potential well of these galaxies makes it easy 
for gas to evaporate or escape in the presence of feedback 
or dynamical heating. We incorporate a simple scheme that 
tracks the gaseous content of galaxies and the effect of gas 
depletion on MBH growth and dynamical evolution. 
(1) 

2.2.1 Halo baryon content 

We track the baryon conte nt of haloes using a method 
similar to that developed by lOkamoto et al] (I2OO8I ). They 
present a simple model that reproduces the results of cos- 
mological hydrodynamical simulations of mass loss from a 
UV background. In this model, haloes form at high red- 
shif ts with the cosm ic mean baryon fraction, = 0.18 
(Du nklev et al.|[2008l ). We define the baryon fraction as the 
ratio between the mass in baryons and the total mass of a 
halo. After reionization (set at 2 = 9 here), low mass haloes 
lose their baryons as a result of th e increased te mperature of 
the intergalactic medium (IGM) (|Gnedinll200ol '). The condi- 
tion for baryon loss depends on the relation of the virial tem- 
perature of the halo, Tvir, to the temperature at which pho- 
toheating and radiative cooling are balanced, Teq. We eval- 
uated Tcq at an overdensity, Aevp = 10^, r epresenting the 
densest, most bound region of the halo (see lOkamoto et al.l 
bOOSl . for a discussion of the model parameters and the al- 
lowed range). If the equilibrium gas temperature in these 
dense regions is greater than the virial temperature of the 
halo, the gas will evaporate out of the halo into the surround- 
ing IGM. In more massive haloes, the dense gas is able to 
cool efficiently, preventing it fro m evaporating. We model 
Tcq using the UV background by iHaardt fc Madaul l|l996l ). 
Halos with Tvir < Tcq(Acvp) have a mass in baryons, Mb, 
which decreases with time, given by: 



My,{t + St) = Mb(t)e . 



(4) 



The evaporation time-scale, t^vp, is given by -Rvir/cs(Aovp), 
where i?vir is the virial radius of the halo, and Cs(Aevp) is the 
sound speed at the evaporation overdensity. This timescale 
corresponds to the time for gas to leave the halo moving at 
the sound speed. 
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We additionally allow for accretion of baryons from the 
IGM onto haloes. If gas at the outskirts of the halo is colder 
than the virial temperature of the halo, the halo accretes 
enough baryons to reach the cosmic mean. Here the tem- 
perature of the accreting gas is approximated as the equi- 
librium temperature of gas with density Pvir/3. This density 
corresponds to the density of gas at the virial radius of the 
halo. The requirement for accretion up to the cosmic mean 
baryon fraction is expressed Tcq(pvir/3) < Tvir- When the 
halo is colder than the gas in the surrounding IGM, no ac- 
cretion occurs. When two haloes merge in our model, the 
resulting halo has a mass in baryons equal to the sum of the 
baryonic masses of the progenitor haloes. 

2.2.2 Massive black hole accretion 

Based on simulations of MBH mergers in gala xies with dif- 
ferent gaseous content (|Callegari et al.ll2009l ). we assume 
that black holes hosted by a baryon rich halo experienc- 
ing a major merger (mass ratio greater than 1:10) accrete 
mass from the host. We here define a halo as baryon rich 
when it has retained more than half of its original baryon 
fraction. Specifically, we choose a baryon fraction thresh- 
old fb > 0.1. In this simple treatment, a merged black 
hole accretes mass according to the central velocity dis- 
persion of the host = 10*(a/200 km s"^)" Mq (e.g., 

iTremaine et al.ll2002l ). We link the central velocity disper- 
sion to the circular velocity of the halo usi ng the empiri- 
cal r elationship log(Vc) = 0.741og((T) -|- 0.8 (jPizzella et al.l 
I2OO5I ). The correlation between central velocity dispersion 
and halo circular velocity has been studied observationally 
for samples of galaxies mostly in the range a > 70kms~^ 
l|Ferraresell2'00a iBaes et al.ll2003l '). Pizzella et al. (2005) ex- 
tend the study to cr ~ 40kms~^, but extending these stud- 
ied to the range of dwarf galaxies is challenging, as the stars 
may not reach the radius where th e halo has its maximum 
circular velocity (e.g., see Fig 5 of iPenarrubia et aLll2008l ). 
Since there is no accepted relationship between Vc and the 
central velocity dispersion that extends down the the dwarf 
galaxy sizes considered here, we use the above relationship 
to determine MBH masses from the halo velocity dispersion. 
For dynamical evolution of the satellites (cfr. Section 2.3.1), 
we use the definition a = Vc/^, as the velocity dispersions 
quoted in the literature that we need to assess our dynami- 
cal modelling are obtained from the 'global' velocity sample 
in each galaxy, not just from some innermost sub-sample. At 
present the only dwarf spheroidal (dSph) satellite with max- 
imum circular velocity constrained by data (and assuming 
constant velocity anisotropy) is Fornax, with g lobal velocity 
disper sion a ~ llkms"^ and Vc ~ 18kms~^ (jWalker et al.l 
l2009al ). consistent with the a = 14 //3 scaling. We note that 
our overall results are robust to changes in the assumed scal- 
ing between Vc and a. We have run models where we also 
adopt the relationship a = Vc//3 to calculate MBH masses, 
and qualitatively all results hold. 

The total mass of the black hole resulting from a com- 
pleted merger is the sum of the merging black hole masses 
and any accreted mass. When a halo with a central black 
hole merges with an empty halo, the central black hole will 
accrete mass after the merger timescale if the halo remains 
baryon rich. A similar scheme has been shown to repro- 
duce observational constraints on MBH evolution (luminos- 
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Figure 1. Top panel: luminosity function of satellite galaxies, 
corrected to match the SDSS DR5 selection. Upper grey curve 
with error-bars: full sample of simulated galaxies surviving tidal 
stripping. Lower red curve with error-bars: galaxies with baryon 
fraction > 7.5x 10~^. Blue circles: observed luminosity function of 
Milky Way satellites. Middle panel: cumulative radial distribution 
of our simulated satellite population (dashed line: from ToUerud 
et al. 2008, fig. 5. 'DR5' distribution). Bottom panel: cumulative 
velocity distribution of our simulated satellite population com- 
pared to the fit derived for satellites in Via Lactea I (lower line, 
Diemand et al. 2007) and Aquarius (upper line, Springel et al. 
2008). Circles: observed dwarf spheroidal satellites of the Milky 
Way. Errorbars are 1— cr Poissonian errors. 



ity function of quasars and Soltan's argument, Mbh — c" 
relationship at z = 0, mass density in MBHs at 2 = 0; 
IVolonteri et al.ll2008l ) for the 'peak3.5' and 'massive seeds' 
models. 



2.3 Dynamical evolution of satellites and massive 
black holes 



As shown by several investigations (e.g. , lMadau et al.ll200j ; 
iTollerud etahlbOOSl : iMaccio et allboiol : fMunoz et al.ll2009[ ) 
the population of dark matter satellites found in numerical 
simulations exceeds by a large degree the number of known 
satellites of the Milky Way, creating a "missing satellite" 
problem. It has been advocated that the solution to this 
problem may lie in a combination of factors. On the one hand 
cooling and star formation are inefficient in the presence of a 
strong photoionizing background, which prevents the devel- 
opment of a conspicuous luminous component (cfr. Section 
2.2.1). On the other hand, tidal stripping of the satellites 
orbiting in the potential of the Milky Way may cause mass 
loss, leading to systems much less massive than they were at 
the time that they merged with their host halo. In Section 
2.3.1 we describe how we model the dynamical evolution of 
satellites in the Milky Way halo in order to derive the prop- 
erties of the sub-population of luminous satellites that have 
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survived until today. In Section 2.3.2 we address the dynam- 
ical evolution of MBH pairs formed during galaxy mergers. 

2. 3. 1 Dynamical evolution of satellites 

After a halo merger, the smaller halo becomes a satellite of 
the more massive system. These satellite haloes evolve in 
the potential well of the host until z = 0, experiencing tidal 
stripping and possibly merging with the host. 

We model the dynamical evolution of satellites 
within the ho st halo potential using analytical techniques 
IVolonteri et al . 2008). For each satellite that merges with 
the main halo of the merger tree, we evolve the satellite-host 
system by integrating the equation of motion of the satellite 
in the gravitat ional potential of th e host (assuming spherical 
NFW profiles; H avarro et al.lll997l ). including the dynamical 
friction term: 



(fr_ GM{r) ^_ A-kG^ InApM^at 



f(x)v 



(5) 



where f{x) = [erf(a;) — {2x/ ^/TT)e~^ ], x = v/V2a, and the 
velocity dispersion a is derived from the Jeans' equation 
for the composite density profile, assuming isotropy (e.g. 
iBinnev fc Tremainel |2008| ). Here M(r) describes the total 
mass of the host within r, p{r) is the total density profile, 
and the second term represents dynamical friction against 
the background. We include the MBH keplerian potential if 
the galaxies host a MBH. The Co ulomb logarithm, In A , in 
equation ((Sjl is taken equal to 2.5 l|Tavlor fc BabulllioOll ). 

The mass of the satellite evolves during the integration 
because of tidal stripping. At every step of the integration 
we compare the mean density of the satellite to the mean 
density of the host halo at the location of the satellite. Tidal 
stripping occurs at the radius within which the mean density 
of the satellite excee ds the density of the galaxy interior 
to its orbital radius l|Tavlor fc Babull l200ll ). We trace the 
orbital evolution and the tidal stripping of all satellites from 
the time when the satellite enters the virial radius of the 
host to z = 0. Satellites that survive until the present time 
provide an analogue of the dwarf galaxy population around 
the Milky Way. 

We compare the circular velocity and radial distribution 
of our satellite p opulation to tha t of V ia Lactea and Aquar- 
ius simulation s (Dicm and et al.l l2007l : iMadau et ahl l2008l : 
ISpringel et al.i .2008) a nd with constraints derived from the 
observed populations |Tollerud et al.l l2008l : IKoposov et al.l 
I2OO8I : IWalsh et al.ll2009^ . Figure [Tl compares our results to 
the literature. We find good agreement with the circular ve- 
locity distribution of satellites in the Aquarius simulation at 
2 = at velocities above Vmax ^ lOkms^^. 

We further derive the luminosity function of the satel- 
lites by assigning a luminosity based on their velocity dis- 
persion. We adopt an emp irical correlations derived by 
iKormendv fc Freeman! (120041 ) for galaxies with velocity dis- 
persion a > 30kms~^: 



log(Ly/LQ) = 7.80 + 51og(cr/30kms"^). 



(6) 



We perform a least square fi t in los(Lv) v s log(a) for low 
a syst ems ( using data fro m Mated (|l998l ). ISimon fc Gehal 
l|2007h and IWalker et all (|2009l )) we find a very sim- 
ilar relationship (log(Lv/LQ) = (7.6 ± 0.6) + (4.3 ± 
0.9) log((T/30kms~^)), albeit with a larger scatter. We 



therefore adopt the fit in Equation |6] for the whole a range. 
To derive the luminosity function we need to set one sin- 
gle parameter, that is the minimum baryon fraction that 
allows star formation. We find the best fit to the observed 
luminosity function by considering as 'luminous' only those 
satellites with a baryon fraction > 7.5 x 10~^ (the accept- 
able range is 6.5 x 10"^ < fb < 8 x 10~^). To comp ute 
the luminosity function, we follow iMadau et al.l (|2008l ) in 
correcting our theoretical sample to match the SDSS DR5 
sample of satellites. First, we correct the sample by a factor 
/dr5 = 0.194 that takes into account the sk y coverage. We 
then apply a joint distan ce-magnitude cut (|Koposov et al.l 
I2OO8I : iTollerud et "aL ''2008) ensuring that we count only satel- 
lites within the SDSS completeness limit: 



47r/D 



1/3 



10 



(-0.6Mv-5.23)/3 



Mpc. 



(7) 



We use the central black holes of this population of 
haloes to provide observational signatures of black hole for- 
mation and growth efficiency. We note that stripping some- 
times (a few per cent of the cases) leaves some 'naked' 
MBHs, as the host loses most of its mass. 



2.3.2 Dynamical evolution of massive black holes 

Along the dynamical evolution of the MBH-|-galaxy sys- 
tems, we further must determine the fate of the MBHs 
they contain (if any). Work by Callegari et al. (2009) has 
shown that the efficiency of black hole pairing is a func- 
tion both of the mass ratio and the baryon fraction of the 
merging haloes. Black holes in merging disc galaxies form 
a bound pair when the mass ratio of the merging galaxies 
is larger than 1:10 and galaxies are gas-rich (i.e., cold gas 
represents 10 per cent of the disc mass). We therefore de- 
fine a major merger to be a merger between haloes with 
mass ratio greater than 1:10. In our models, minor mergers 
(mass ratio less than 1:10) do not lead to efficient black 
hole pairing and mergers. When a major merger occurs, 
the final fate of the MBHs depends on the baryon content 
of the host. There is insufficient dynamical friction in gas- 
poor galaxies to lead to efficient MBH pairing. We assume 
that black hole mergers stall when occurring in a baryon 
poor halo (/{, < 0.1). Unless stalled, we assume that black 
holes merge within the merger timescale of their host galax- 
ies. For major mergers, this timescale is well approximated 

by Tinerge /rdyn ^ 0.4(Mho.t /M^at/ / ln( 1 + Mhost/M^at), aS 

shown bv lBovl an-Kolchin et al. (2008), where Tdyn is the dy- 
namical time at the virial radius. This times cale represents a 
lower limit to the MBH-MBH merger time (jBegelman et al.l 
Il980l ). although the assumption that MBHs merge within 
the merger timescale of their hosts is li kely for MBH pairs 
formed after gas rich galaxy mergers (|Escala et all l2005l : 
iDotti et al]|2007D . 



3 BLACK HOLE OCCUPATION FRACTION 
AND MASS SCALING 

We study the black hole population of satellites at z = 
to find signatures of the black hole properties and seed effi- 
ciency. We show the black hole occupation fraction (BHOF) 
of the satellite population at z = in Figure (2] Massive 
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satellites today were likely to have massive progenitors at 
high redshifts, meaning that they are more likely to host 
central black holes, given that a minimum mass threshold is 
required in all our models (mimicking the necessity of a deep 
potential well). The higher the threshold mass, the more the 
BHOF is expected to decrease with the mass (or velocity dis- 
persion) of today's host. As expected, the more biased sce- 
nario 'peaks. 5' leads to a lower overall BHOF than 'peakS'. 
The 'massive seed' scheme produces fewer black holes than 
either 'peakS' or 'peak3.5'. This is because of the stricter 
conditions needed for MBH formation: seeds form only in 
very massive halos (Afh — 10* M0, compared to the much 
lower mass threshold for the Population III remnant case, 
Mh ~ 10'' M© at 2 = 15 for 'peakS' and Mh ~ 7 x 10^ M© 
for 'peaks. 5') with low angular momentum (low spin param- 
eter). Given that spin parameters appear to be distributed 
lognormally (with mean Aspin = 0.05 and standard deviation 
a\ = 0.5) in halos extracted from cosmological simulations, 
typically only ~10 per cent of the haloes in the allowed mass 
range can form a central seed MBH, for the parameter choice 
discussed in section 2.1. The results shown here are largely 
insensitive to the efficiency of black hole merging and ac- 
cretion. The merger and accretion efficiency primarily affect 
the mass of the central black holes, not the presence or ab- 
sence of one. The present-day BHOF is therefore a sensitive 
probe of the efficiency of black hole formation in haloes at 
high redshifts. Observations of satellites at relatively high 
velocity dispersions might therefore distinguish between the 
massive seed and Population III models presented here, al- 
though this prospect depends on the ability to detect MBHs 
where they exist (see Section 4). 

Unlike the BHOF, the masses of MBHs within galaxies 
are in general sensitive to the merger efficiency and accre- 
tion. However, in the environment we are investigating we 
expect MBH growth to be largely inefficient. In our model, 
black holes accrete mass when the host halo experiences 
a major merger and remains baryon rich over the merger 
timescale. Gas-rich major mergers are rare for the progeni- 
tors of our satellite galaxies, so MBH growth through merg- 
ers and accretion is suppressed. 

In a baryon poor halo, MBH-MBH mergers will not 
complete, accretion does not occur, and the black hole mass 
does not change. Even before reionization, when all galaxies 
are gas-rich, two factors contribute to the limited growth of 
MBHs. On the one hand, the rarity of seeds causes MBHs to 
evolve mostly in isolation. On the other hand, MBH growth 
through accretion at high redshifts is negligible in our model. 
Accreting MBHs grow according to the velocity dispersion 
of the host halo, but the seed masses we consider exceed 
the MBH mass expected for small haloes. For satellites with 
velocity dispersions similar to those of the Milky Way dwarf 
galaxies, Population HI remnant black holes grow to no more 
than about an order of magnitude larger than the initial 
seed mass at z = 0. Massive seed masses remain effectively 
unchanged. This is an important result, independent of the 
formation scenario. MBHs in dwarf galaxies indeed provide 
a clear indication of the initial seed properties. 

This can be seen in Figure [2l where we show the ex- 
pected relationship between MBHs and their 'dwarf hosts. 
At (J < 20 — 40kms~^ there is no correlation between 
MBH masses and velocity dispersion. This is due to the 
very limited mass growth of MBHs hosted in low-mass satel- 
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Figure 2. Top: fraction of galaxies, at a given velocity disper- 
sion, which host a central MBH (occupation fraction, BHOF). 
Long-dashed green histogram: 'peak3.5' case. Short-dashed green 
histogram: 'peak3' case. Thick lines: all sateUites. Thin lines: lu- 
minous satellites. Sohd red histogram: 'massive seeds'. Bottom: 
the Mbh— velocity dispersion (cr) relationship for MBHs in satel- 
lites. We here show the results of a suite of ten realisations. Stars: 
'massive seeds'. Circles: 'peak3.5'. Triangles: 'peakS'. 



lites. These satellites have a very quiescent merger history 
(very few satellites experience any major merger at all at 
z < &), causing their central MBHs to starve and remain 
near the original seed mass, creating an asymptot ic hori- 
zontal 'plume' (see also I Volonteri fc Nataraiar]l2009l i. Addi- 
tionally, a secondary role is played by mass stripping, that 
depends on both the mass of the host and the satellite, but 
also on the orbital parameters - stripping being much more 
effective for radial orbits where satellites plunge deep into 
the potential well of the host, where densities are higher. A 
severely stripped system has a much larger MBH-to-host ra- 
tio. The final MBH-to-host relationship is therefore a combi- 
nation of nature (formation mechanism) and nurture (MBH 
feeding and dynamical evolution of the host). 



4 OBSERVATIONAL PROSPECTS 
4.1 Dynamical measurements 

The presence of an MBH can be tested dynamically if the 
region where the gravitational potential of the black hole 
dominates the gravitational potential of the host can be 
resolved. This region is referred to as the sphere of influ- 
ence of the black hole. We adopt here the definition of the 
sphere of influence as the region within which the enclosed 
mass in stars and dark matter equals twice the MBH mass. 
The radius of the sphere of influence is therefore defined as: 
M{r < Anf) = 2 X Mbh- 

The lower panel in Figure [3] plots stellar velocity dis- 
persion against Rinf, estimated for massive seed and Pop 
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Figure 3. Top panel: number of stars within a given projected 
radius in Fornax. Lower curve; all Fornax member stars for which 
velocities are currently available in the published kinematic sam- 
ples of Walker et al. (2009a). Upper curve: all Fornax target can- 
didates, including unobserved stars, that are sufficiently bright 
{V < 21.5) for velocity measurements with existing 6m - 10m 
telescopes. Bottom panel: relationship between velocity disper- 
sion ((t) and radius of the sphere of influence of MBHs — defined as 
the sphere that encompasses 2 X Mbh - for ten halo realisations. 
Stars: 'massive seeds'. Circles: 'peak3.5'. Black triangles: Milky 
Way satellites: we assume that the MBH sits on the Mbh ~ 
relationship, and the density profile is cored. Grey triangles: as 
above, but for an NFW profile. Black dots: we assume a fixed 
MBH mass, 10^ Mq and a cored profile. Grey dots: as above, but 
for an NFW profile. 



also iBattaglia et al.ll200^ ). and 2) all Fornax target candi- 
dates, including unobserved stars, that are sufficiently bright 
(V^ < 21.5) for velocity measurements with existing 6m - 
10m telescopes. The latter profile represents the largest sam- 
ples that are possible at present; unfortunately, these would 
include fewer than 5 stars within the spheres of influence 
estimated for the classical dSphs. Thus even for the bright- 
est nearby dSphs, the detection of any MBH must await the 
next generation of 20-30m telescopes, which may increase 
kinematic sample sizes by more than an order of magnitude. 

Finally, we consider the spheres of influence due to MBH 
masses of ~ 10^ Mq , a mass of order of the upper lim- 
its derived for the 'massive se ed' scenario (we reca l l here 
that the formalism proposed bv lLodato Sz NataraianI (|2006l ) 
and adopted here gives the total mass available for forming 
an MBH seed. Most like ly the final seed mass is smaller, 
see iBegelman et al.ll200^ ). For such masses, the expected 
spheres of influence reach < 50 pc for the observed dsphs. 
For Fornax, the expected value of 7?inf ~ 30 pc encloses 10 
stars in the existing velocity sample, and 25 stars in the list 
of current target candidates. If all these stars are observed, 
the resulting sample may help diagnose whether Fornax has 
an MBH of mass ~ lO'^ Mq. 

Figure [3] suggests that MBHs generated as 'massive 
seeds', having larger masses and a larger J?inf, would be 
the most favourable scenario from an observational point 
of view. However, their typical BHOF is lower, being always 
below 40 per cent and decreasing to less than a per cent for 
'true' dwarf galaxy sizes (a ~ 5-15 km s~^; I Walker et al.l 
l2009bl ). From Figure [1 the typical BHOF for a massive seed 
MBH in a dwarf galaxy like Fornax is ~ 2 per cent, mean- 
ing that ~ 50 Fornax-like dwarfs would need to be observed 
to have significant probability of finding an MBH. Popula- 
tion HI remnants have instead a higher BHOF, but their 
masses have not grown much since formation, making their 
detection harder. 



Ill cases. Overplotted are values estimated for the eight 
'classical' dwarf spheroidal satellites of the Milky Way, for 
which line-of-sight velocities have b een measured for u p 
to a few thousand stars per gala:xy () Walker et al.l lioOQal ) . 
For these object s we adopt the stella r velocity dispersion 
measurements of I Walker et al.l (|2009bl ). and then adopt an 
MBH mass from the mass-velocity dispersion relationship 
iTremaine et al.ll2002l ). In order to calculate the sphere of in- 
fiuence for the real dSphs, we consider the best-fitting cored 
and cuspe d (NFW) mass profil es with parameters listed in 
Table 3 of IWalker et"all (|2009bl ). If we assume the observed 
dSphs have NFW dark matter profiles, then according to the 
M-sigma relation, their MBH masses correspond to spheres 
of influence of 1 < 7?inf < 2 pc. If the dSph dark matter 
profiles are cored, then the spheres of influence are a few 
times larger, 3 < Rinf < 7 pc. 

In order to evaluate prospects for detecting kinematic 
signatures from such MBHs in real dSphs, the upper panel 
of Figure O indicates the number of spectroscopic target 
stars within a given projected radius in Fornax, the most 
luminous dSph satellite of the Milky Way. Curves indicate 
the cumulative surface brightness profiles of 1) all Fornax 
member stars for which velocities are currently available in 
the published kinematic samples of Walker et al. 2009a (see 



4.2 Gravitational waves 

The detection of gravitational waves from a MBH in a dwarf 
galaxy undergoing a merger with another black hole is an- 
other possible probe. The Einstein Telescope, a proposed 
third-generation ground-based gravitational-wave (GW) de- 
tector will be able to detect GWs in a f requency range reach- 
ing down to ~ 1 Hz (|Freise et al.ll2009l ). Sources with masses 
of hundreds or a few thousand solar masses will generate 
GWs in this frequency range, which is out of reach of LISA 
or the current ground-based detectors. Since dwarf galax- 
ies have a very quiet merger history, we would not expect 
many MBH-MBH mergers involving dwarf galaxies at the 
present epoch. However, gravitational waves may also be 
generated in dwarf galaxies by mergers between the cen- 
tral MBH and stellar remnants (stellar mass black holes) 
in the centre of the dwarf. It is estimated that the Ein- 
stein Telescope might detect as many as one thousand of 
these inspiral events in globular clusters, assuming a rela- 
tively high intermediate mass black hole occupation frac- 
tion in the clusters (Gair et al. , 2009i). The same mechanism 
would operate in dwarfs, and the predicted rate of mergers 
from this channel scal es with the core stellar density, n, as 
n'*''^ (|Gair et al.ll2009l ). The core stellar densities in nearby 
dwarf galaxies are typically very low, however, e.g., the esti- 
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mate for Fornax is ~ 10~^pc~^ (|M atecilll99^) and for Sagit- 
tarius is ~ 10~^pc~^ (jMaiewski e t al. 20051). These densities 
are several orders of magnitude lower than the typical core 
densities in globular clusters, which are ~ lO^'^pc"^. 

We estimate the number density of MBHs in dwarfs 
via recent theoretical works that study the population of 
dwarfs as sate l lites o f the Milky Wav jReed et af] l2005l : 
iDiemand et~al] l2007l : ISpringel et ahl bOOSl 'l. These simula- 
tions suggest that the number of satellites per halo has the 
following form: 



iV(> Vs: 



( Vs. 



\ Whost 



(8) 



where i^sat and Uhost are the majdmum circ ular velocity of 
the sa tellite and the host halo, respectively. ISpringel et al.l 
(|2008ll find iV. = 0.052 and a = -3.15. If we extrapolate 
the Mbh — o" correlation to 10^-10^ Mq BHs and assume 
an isothermal galaxy profile, then Waat ~ 10-20 km s~^. 
With this formalism we obtain the number of satellites in the 
interesting mass range per dark matter halo (Nsa^t), where 
the mass of the halo is uniquely determined by its maximum 
circular velocity. The number density of dark matter halos 
can be easily obtai ned by integrating the m odified Press & 
Schechter function (|Sheth fc Tormen Il999l ) which provides 
the mass function of halos, dn/dMh- Therefore we estimate 
a number density of satellites (per comoving cubic Mpc) as: 



dn 
dMh 



Ns^tiMh)dMh ~ IMpc" 



(9) 



adopt ing the normalization proposed by ISpringel et al.l 
l|2008t ). We now have to correct for the fact that not all 
dwarf galaxies are likely to host an MBH. Based on the 
models presented here, a fraction ~ 0.01-0.1 of dwarfs host 
an IMBH with mass ~ 10^-10^ Mq. Furthermore, we must 
accou nt for the fact that man y of these satellites do not form 
st ars teovill fc Ricotti|[2009l . and references therein). Based 
on accio et al.l()2010[ ). we estimate that a fraction /* ~ 0.8 
of dwarfs with masses > 2 x 10*^ Mq formed stars and will 
contain stellar mass BHs that can merge with the central 
IMBH. The final estimate for the number density of dwarfs 
hosting an MBH that could be Einstein Telescope sources is 
~ 0.008-0.08 Mpc""^. This number density is approximately 
an order of magnitude small er than the ~ 0. 3 Mpc""^ as- 
sumed for globular clusters in iGair et all (|2009l ). Combining 
the core stellar density and the number density scalings, we 
conclude that the rate of GW detections from this channel 
is likely to be <C 1 per year. Therefore, although it is not 
inconceivable that third generation GW detectors will de- 
tect events originating in dwarf galaxies, any events will be 
serendipitous. 

An alternative way of detecting MBHs lurking in dwarf 
galaxies would be via their emission when accreting sur- 
rounding material, either from a companion star or gas avail- 
able as recycled material via mass loss of evolved stars (Dotti 
et al. in preparation). 



5 CONCLUSIONS 

We study the black hole population of dwarf galaxies in a 
Milky Way-type halo to find observational signatures of the 
MBH seeding and growth processes. Dwarf galaxies have 



a quiet merger history, leading to an MBH population at 
z = that has not evolved significantly from the seed pop- 
ulation at high redshift. We consider two MBH formation 
mechanisms: a 'massive seeds' model where seed form via 
gas-dynamical collapse and a Population HI remnant model. 
We seed the progenitor halos of a Milky Way-type halo with 
MBH seeds and follow their evolution until z — 0. Incom- 
plete halo mergers at low redshifts produce a satellite popu- 
lation consistent with observated and simulated results. We 
study the MBH population of these satellites and find: 

• Most dwarf galaxy MBHs have not grown significantly 
since seed formation at large redshift. Dwarf galaxies have 
low masses and tend to be baryon-poor, suppressing MBH 
growth through mergers and accretion. Measurements of 
MBH masses in dwarf galaxies would directly probe the seed 
masses. 

• The 'massive seeds' model produces rare but large seeds 
(M ~ 10'=' Mq). The resulting MBH population at z = oc- 
cupies less than a few per cent of typical dwarf galaxies. 
These MBHs are sufficiently massive that they may be de- 
tectable from the central stellar kinematics in bright dwarfs 
such as Fornax. 

• The Pop HI remnant model produces many small seeds 
(M = 10'^ Mq). This MBH population is more abudant at 
low redshifts than in the 'massive seeds' model, but the low 
mass MBHs cannot currently be detected. The detection of 
this population must await the next generation of telescopes. 
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